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Abstract 

We consider the motion of a spinning relativistic particle in external electromagnetic 
and gravitational fields, to first order in the external field, but to an arbitrary order in 
spin. The correct account for the spin influence on the particle trajectory is obtained 
with the noncovariant description of spin. Concrete calculations are performed up to 
second order in spin included. A simple derivation is presented for the gravitational 
spin-orbit and spin-spin interactions of a relativistic particle. We discuss the gravimag- 
netic moment (GM), a specific spin effect in general relativity. It is demonstrated that 
for the Kerr black hole the gravimagnetic ratio, i.e., the coefficient at the GM, equals to 
unity (as well as for the charged Kerr hole the gyromagnetic ratio equals to two). The 
equations of motion obtained for relativistic spinning particle in external gravitational 
field differ essentially from the Papapetrou equations. 
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1 Introduction 



The problem of motion of a particle with internal angular momentum (spin) in an external 
field consists of two parts: the description of the spin precession and the account for the 
spin influence on the trajectory of motion. To lowest nonvanishing order in c~ 2 the complete 
solution for the case of an external electromagnetic field was given more than 70 years ago |3J . 
The gyroscope precession in a centrally symmetric gravitational field had been considered 
to the same approximation even earlier ||. Then, essentially later the spin precession was 
investigated in the case of the gravitational spin-spin interaction M. The fully relativistic 
problem of the spin precession in an external electromagnetic field was also solved more than 
70 years ago || , and then in a more convenient formalism, using the covariant vector of spin, 
inRef. [|. 

The situation with the second part of the problem, which refers to the spin influence on the 
trajectory, is different. Covariant equations of motion for a relativistic spinning particle in an 
electromagnetic field were written in the same paper [|J, and for the case of a gravitational 
field in Ref. ||. These equations have been discussed repeatedly from various points of 
view in numerous papers [7-15]. The problem of the spin influence on the trajectory of a 
particle in external fields is not of a purely theoretical interest only. It attracts attention 
being related to the description of the motion of ultrarelativistic particles in accelerators [[Uj 
(see also recent review |17|j ). Then, there are macroscopic objects for which internal rotation 
influences their motion in an external gravitational field. We mean Kerr black holes. This 
problem is of importance in particular for the calculation of the gravitational radiation of 
binary stars. In this connection it was considered in Refs. [18-21]. However, when turning 
to these calculations, we found out [22 that the equations of motion with the account for 
spin to the lowest nonvanishing order in c~ 2 , used in these papers, even in the simpler case 
of an external field lead to results which differ from the well-known gravitational spin-orbit 
interaction. The problem is essentially related to the correct definition of the center-of-mass 
coordinate. It turned out moreover that the Papapetrou equations as well do not reproduce 
in the same c -2 approximation the result for the gravitational spin-orbit interaction ascending 
to Ref. 0. This discrepancy was pointed out long ago in Ref. |23[], however its explanation 



suggested there does not look satisfactory (see Ref. f22j). 

In the present work we obtain equations of motion of a relativistic particle with the 
noncovariant description of spin. They agree with well-known limiting cases. Though for 
external electromagnetic field such equations in the linear in spin approximation have been 
obtained previously |U| (see also Ref. [0), we would like to start with comments related to 



this approximation in electrodynamics. 
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2 Covariant and noncovariant equations of motion 
of spinning particle in electromagnetic field 

The interaction of spin with external electromagnetic field is described, up to terms on the 
order of c~ 2 included, by the well-known Hamiltonian 

ff=-^sB + ^ S [pxE]. (1) 
2m 2m 2 [P J K J 

Let us emphasize that the structure of the second, Thomas, term in this expression not only 
is firmly established theoretically, but has been also confirmed with high accuracy experi- 
mentally, at any rate in atomic physics. The force corresponding to the Hamiltonian ([!]) 
equals 

f ™= £ sB - +£ ^(si Exs i-- s i vxE -]) ; < 2) 

here and below coma with index denotes a partial derivative. 

Let us try to construct a covariant equation of motion accounting for spin, which would 
reproduce in the same approximation the force (Q). A covariant correction to the Lorentz 
force eF^ u u ul linear in the tensor of spin S^ u and in the gradient of the tensor of electro- 
magnetic field -F^a , may depend also on the 4- velocity u^. Since u^u^ = 1, this correction 
should satisfy the condition = 0. From the mentioned tensors one can construct only 
two independent structures which satisfy the last condition. The first one, 

ri^F vX>K S vX - F Xu , K u*S x »u^ 

reduces in the approximation of interest to 

2s(B, m - [v x E, m ]), 

and the second, 

u x F Xv ^S»^ 

reduces to 

|[sxE] m . 

(Let us note that the structures with the product F UKjX S kX reduce to the two presented 
expressions, due to the Maxwell equations and the antisymmetry of S kX .) Obviously, no 
linear combination of the two presented structures can reproduce the correct expression (|2]) 
for the spin-dependent force. In a somewhat less general way it was demonstrated in Ref . |22| . 
It was pointed out in the same paper, that the coordinate entering the covariant equation does 
not coincide with the usual one. Therefore, to obtain the correct c~ 2 approximation to the 
covariant equation of motion one needs to perform an additional redefinition of coordinate: 

r — > r v x s. (3) 

2m 
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In the case of spin 1/2 this redefinition is closely related to the Foldy-Wouthuysen transfor- 
mation pi| . The generalization of this substitution to the case of arbitrary velocities was 
recently suggested in Ref. JT7 . 



Meanwhile, the correct equations of motion in electromagnetic field including spin to first 
order are known for sufficiently long time |T!| . Let us recall that the initial physical definition 
of spin refers to the proper frame of the particle. This is the 3-dimensional vector s (or 3- 
dimensional antisymmetric tensor) of the internal angular momentum defined in this frame. 
The covariant vector of spin (or the covariant antisymmetric tensor S^ u ) are obtained 
from them just by the Lorentz transformation. By the way, in this approach the constraints 
u^Su = 0, u^Snu = are valid identically. The precession frequency for spin s at an arbitrary 



velocity is well-known (see, for instance, Ref. |25|j): 
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7 + 1 
1 



v(vB) - v x E 



v x E 



7 



(4) 



where 7=1/ Vl — v 2 . Naturally, the corresponding interaction Lagrangian (the Lagrangian 
description is here somewhat more convenient than the Hamiltonian one) equals 



2m 



+2 



(9-2) 



B 




v(vB) - v x E 



_7 7 + 1 

The force, we are looking for, is given by the usual relation 

d 



(V 



dt 



:VvL, 



(5) 



(6) 



and the equation of motion for spin in general form is 

s = -{L,s}, (7) 
where {... , ...} denotes the Poisson brackets, or 

s = -i[L,a] (8) 

in a quantum problem. 

In the conclusion of this section, let us discuss the following question. In fact, it is far 
from being obvious how meaningful are at all the discussed spin corrections to the equations 
of motion of elementary particles, say, electron or proton. According to the well-known 
argument by Bohr (see Ref. f26|), an additional Lorentz force due to the finite size of the 
wave packet of a charged particle and to the uncertainty relation, exceeds the corresponding 
component of the Stern-Gerlach force. Meanwhile, it was proposed long ago to separate by 
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polarizations a beam of charged particles in a storage ring through the spin interaction with 
external fields [p7 



Though this proposal is being discussed rather actively (see review 
it is not clear up to now whether it is feasible in principle. In this connection we would like 
to note that, as distinct from the Lorentz force e(E + v x B), the component of the spin force 
(0) directed along the velocity, grows with energy, its asymptotic behaviour being 



7V (I +VV ) ^ S 



o-2 

B - v x E + v(vB) 



This growth of the spin force with energy gives at least in principle the possibility to circum- 
vent the limitation due to the uncertainty relation. 



3 Equations of motion of spinning particle 
in electromagnetic field. General formalism. 
Effects linear in spin 

In this section we point out a general approach to the derivation of the equations of motion 
in external electromagnetic field to an arbitrary order in spin. We will reproduce here in 
passing the known result (Q). 

We derive the Lagrangian of the spin interaction with an external field from the scattering 
amplitude 

- e JM„ (9) 

of a particle with spin s off a vector potential A^. Due to the arguments presented in the 
end of the previous section, the account for the effects nonlinear in spin, which we are first 
of all interested in, may have physical meaning only in the classical limit s ^> 1. Just this 
approximation will be mainly used below. 

The matrix element J M of the electromagnetic current operator between states with mo- 
menta k and k! can be written (under P and T invariance) as follows (see, for instance, Refs. 

^ = ^ {p^F e + Z^q" F m } ip(k). (10) 

Here p M = (k' + fc) M , g M = {k 1 - 

The wave function of a particle with an arbitrary spin ip can be written (see, for instance, 



Ref. J25[, §31) as 

nh = 



♦4 . E l mi 



Both spinors, 
and 



r <r Q i a2 •• a p \ 

^$102- $„ 
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are symmetric in the dotted and undotted indices separately, and 



p + q = 2s. 

For a particle of half-integer spin one can choose 

1 1 

p = s+ 2' q = S ~ 2 ' 

In the case of integer spin it is convenient to use 

p = q = s. 

The spinors £ and r] are chosen in such a way that under reflection they go over into 
each other (up to a phase). At p ^ q they are different objects which belong to different 
representations of the Lorentz group. If p = q, these two spinors coincide. Nevertheless, we 
will use the same expression ([TT]) for the wave function of any spin, i.e., we will also introduce 
formally the object rj for an integer spin, keeping in mind that it is expressed in terms of 
£ . This will allow us to perform calculations in the same way for the integer and half-integer 
spins. 

In the rest frame both £ and rj coincide with a nonrelativistic spinor £o> which is symmetric 
in all indices; in this frame there is no difference between dotted and undotted indices. The 
spinors £ and f] are obtained from £o through the Lorentz transformation: 

£ = exp{£0/2K o ; rj = exp{-E0/2K o . (f 2) 

Here the vector </> is directed along the velocity, tanh0 = v; 



v p+q 

s = 5>- £ 



i -i 



i=l i=p+l 

and crj acts on the ith index of the spinor £ as follows: 

C-i£o = ((TijaiPi (£o)....&... • (13) 

In the Lorentz transformation ( |T2"D for £, after the action of the operator S on £o the first p 
indices are identified with the upper undotted indices and the next q indices are identified 
with the lower dotted indices. The inverse situation takes place for 77. 

Let us note that in an external field the components of velocity v (and together with them 
the components of </>) do not commute, generally speaking. However, to the approximation 
linear in external field, we are interested in, one can neglect this noncommutativity which is 
itself proportional to the field. Moreover, we are mainly interested in the classical limit of 
the final result where such commutators are inessential being proportional to an extra power 
of H. Therefore, v and (f) will be treated as usual numerical parameters. 

Then, 

i> = ^ f 7o = ^ ( j J ; 
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here / is the sum of unit 2x2 matrices acting on all indices of the spinors £ and rj . The 
components of the matrix £^ = — "E Ufl are: 



( -Sn 

JOn ~[ o s n 



(14) 



The scalar operators F e>m depend on two invariants, t = q 2 and r = (S^q^) 2 . The 
covariant vector of spin S M is defined, e.g., for the state with momentum and is obtained 
through the Lorentz transformation from the spin vector (0, s) in the rest frame: 

S" = (5b,S), S = ^, S = s+ k( ^ s) . . (16) 

m m[e + m) 

In the expansion in the electric multipoles 



Me 



F e {t,T) = Y,f e , 2n {t)T n 



n=0 



the highest power N e equals obviously to s and s — 1/2 for integer and half-integer spin, 
respectively. In the magnetic multipole expansion 



Nr, 



Fm(t,T) = Y, fm,2n(t)r n 
n=0 

the highest power N m constitutes s — 1 and s — 1/2 for integer and half- integer spin. It can 
be easily seen that 

/e, (0) = 1, / m , (0) = |. 

Let us note at last that we have chosen the noncovariant normalization for the amplitude 
being interested in the Lagrangian referring to the world time t, but not to the proper 
time t. 

Let us reproduce now in this approach the well-known result (^) for the case of a constant 
external field. We will start with the terms proportional to g-factor. The corresponding term 
in the scattering amplitude can be written as 

e £ $ {[exp{S0/2}(sB) exp{-£0/2} + exp{-S0/2}(sB) exp{£</>/2}] 
+ i [exp{£</>/2}(£E) exp{-£0/2} - exp{-£</>/2}(£E) exp{£0/2}]| f . (17) 
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It is essential that in the case considered, that of the constant external field, one may put 
k' = k, v' = v, </>' = </>, since q = k' — k corresponds to the field gradient. 
The further calculations use the known identity 



exp{A}Bexp{-A} = B + - [A, B] + - [A, [A, B] 



+ ... 



as well as the relations 

[Sj, Hj] = 4i€ij k s k , [Sj, Sj] = ieij k H k ; (18) 

cosh0 = 7, sinh0 = t>7. (19) 
After sufficiently simple algebraic transformations expression (|T7j) reduces to 



eg 

— s 

2m 



B — v(vB) - v x E 



7 + 1 

Let us go over now to the contribution of the convection term 



(20) 



- ^(k'Mk)^. (21) 
We rewrite the product of exponents in the expression 

$(k')i/>(k) = ^o t [exp{S072}exp{-S0/2} + exp{-S072}exp{S0/2}]e o (22) 

as 

exp{E0'/2} exp{-£</>/2} 

= J]exp{cr072}exp{-o-0/2} J]exp{-tr072}exp{<r0/2}. (23) 
p g 

Let us consider a typical factor in this formula: 

exp{<7072} exp{-<7</>/2} = cosh(<///2) cosh(0/2) - (n'n) whh{<// /2) sinh(0/2) 

+ er [n' sinh(072) cosh(0/2) - n cosh(<///2) sinh(0/2)] (24) 

- 1 (cr[n' x n]) sinh(072) sinh(0/2); 

here n' = v'/v', n = v/v. Since we are interested in gradients only as long as they enter 
together with spin, in the first term cosh(072) cosh(0/2) — (n'n) sinh(<///2) sinh(0/2) we 
put 0' = 0/2, n' = n, after which this term turns to unity. Then, we are discussing the 
interaction linear in spin, so that the product (|23| ) reduces to 

1 + 53 [n'sinh(072) cosh(0/2) - ncosh(<///2) sinh(0/2)] - 2z (s[n' x n]) sinh(072) sinh(0/2). 



When substituted into formula fl22|), the terms proportional to 53 cancel. Restricting then to 
the terms linear in q, we reduce the spin-dependent part of (pi]) to 

i(s[k x q]) 
2e m(e + m) 
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Let us note further that since p^q^ = , the identity takes place 

TfqaAy, = p^iqaA^ - q^A a ) = p»iF aii . 



(25) 



Now we can put p^ — * 2mu tl , where is the 4-velocity. In result we arrive at the following 
expression: 



2m 



2 1 



7, 



B — v(vB) — vxE . (26) 

7+1 7+1 J 

Expressions fl2"0|), ( |2"6j ) give in sum formula (|5[). Thus, we have reproduced the known result 
for the interaction linear in spin, starting from the relativistic wave equation for an arbitrary 
spin. 

Below we will use repeatedly identities of the type fl25|). In the classical language such 
a transformation corresponds to omitting in a Lagrangian (or adding to it) a total time 
derivative. Indeed, 



27 



7 



d_ 

dt 



+ vV 



7 



d_ 
dt' 



4 Equations of motion of spinning particle 

in electromagnetic field. Effects quadratic in spin 

Let us go over now to the interaction of second order in spin. The "bare", explicit quadrupole 
interaction present in expressions (0), fllCf) is 



e P -f e , 2 {S a q a ) 2 A,. 



(27) 



Using identity ( p5| ) and relations (0), we reduce it to 



- e f e ,2 7 { (vs) 



d 7 
dt + 7 + 1 



VV) 



+ i(.V) 



fsE) 



7 



7 



(sv)(vE) + (s[v x B]) 



Neglecting the total time derivative d/dt + vV, we arrive at expression 



'2s 



~ e fe,2 



(sV) - 



7 



7 + 1 



(vs)(vV) 



(sE) - 



7 



7 + 1 



(sv)(vE) + (s[v x B]) 



(28) 



Using the Maxwell equations and adding a total derivative in t, one can demonstrate that this 
expression has such a structure that the tensor SiSj in it can be rewritten in the irreducible 
form: SiSj —>■ SiSj — (l/3)5jjS 2 . Now it is clear from the nonrelativistic limit of formula (|28| ) 
that it describes indeed the interaction with external field of the quadrupole moment 



Q 



-2ef e , 2 s(2s-l) 



(29) 



In the asymptotics, as 7 — > 00, the interaction (|28| ) tends to a constant 

Us = -e/ e ,2 [(sV) - (vs)(vV)] [(sE) - (sv)(vE) + (s[v x B])]. (30) 

It is well-known that even in the absence of the bare quadrupole term, i.e., at f e ^ = 0, a 
quadrupole interaction arises in the nonrelativistic limit due to the convection and magnetic 
terms in interaction (|S|). The value of this, induced quadrupole moment at an arbitrary spin 
of the particle was obtained in Ref. |29| 0: 

0l - -e(a-l) f— V i *' integer spin, 

6[9 ' \mc) {8-1/2, half-integer spin. 

We have singled out explicitly in this formula the Planck constant % to demonstrate that the 
induced quadrupole moment Q± vanishes in the classical limit H — > 0, s — > 00, hs — > const. 
Therefore, the interaction of second order in spin proportional to Q\ does not influence in 
fact equations of motion of a classical particle (though it plays a role in atomic spectroscopy 

II)- 

Still, the convection and magnetic terms in expression (|9|) induce an interaction of second 
order in spin which has a classical limit and is therefore of interest for our problem. It is 
convenient here to start with the convection current interaction. Let us come back to formula 
(HI). Again we put in it 

cosh(0'/2) cosh(0/2) - (n'n) sinh(0'/2) sinh(0/2) = 1. 

And in other terms, linear in cr, we keep the first power of q — > —ih'V only, in the hope 
that in the final answer, in the product (^), h will enter in the combination hs — ► const. 
Nevertheless, these terms by themselves are small as compared to unity, so that in the classical 
limit expression ( Ell ) can be rewritten as 

exp {<r [n'sinh(072) cosh(0/2) - ncosh(0'/2) sinh(0/2)] - % (tr[ri x n]) sinh 2 (0/2)} . 

It can be easily seen that in the product (|23| ) the operators cr, entering with 

n' sinh(0'/2) cosh(0/2) - ncosh(0'/2) sinh(0/2) , 

combine in the resulting exponent into the operator XI which vanishes in the classical limit. 
In this limit only those operators cr survive which enter with [n' x n] sinh (0/2); they 
combine into 2s. Thus, in the classical limit the product (|23| ) reduces, with the account for 
the second identity (0) to 

ex p{^I («[vx V])|. (32) 

Let us note that that the action of the operator fl32p on any function of coordinates, 
whether it is a vector potential or field strength, reduces to the shift of its argument: 

1 7 

r^rH s x v . 

m 7 + 1 

3 The authors of this paper overlooked a misprint in the answer for the induced quadrupole moment, in 
result its value presented in Ref. (29) is two times smaller than the correct one 
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Curiously, just this substitution is suggested in Ref. [|T7[] for the transition from the covariant 
equations linear in spin to noncovariant ones. Its particular case in the c -2 approximation is 
formula (H). 

Now, taking into account the second term in the expansion of the exponential function 
fl32|) and using again the identity fl25p, we easily obtain the following expression for the 
quadratic in spin interaction arising from the convection current: 



7 



2m 2 7 + 1 



(s[vx V]) 



ir 7 + 1 



7 



sv)(vB) - — ^(sfvxE]) 



(33) 



Let us go over now to the contribution into the discussed effect due to the magnetic 
moment. The term in the scattering amplitude due to the magnetic moment, we are interested 
in, can be conveniently written now as 

^eo t {[exp{S072}exp{-S0/2}exp{S0/2}(sB)exp{-S0/2} 
+ exp{-S0'/2} exp{£0/2} exp{-£</>/2}(sB) exp{£</>/2}] 



+ - [exp{£072}exp{-£(/>/2}exp{£0/2}(£E)exp{-£0/2} 
- exp{-£0'/2} exp{£(/>/2} exp{-£0/2}(£E) exp{£0/2}]} £ - 



(34) 



Using in this case the first term in the expansion of the exponential function (0), we arrive 
at the following expression for the contribution proportional to the magnetic moment: 



eg 7 



2m 2 7 + 



- (s[vx V]) 



fsB) 



7 



7 + 1 



(sv)(vB) 



S V 



xE]) 



(35) 



The total result for the induced interaction, quadratic in spin, is 



7 



J 2s 



2m 2 7 + 1 



(s[vx V]) 



9-1 + 



7, 



fsB) 



(9 



7 



7 + 1 



(sv)(vB) 



9 



7 



7 + 1 



(s[vxE]) 



(36) 



Let us note that in the nonrelativistic limit the induced interaction with magnetic field 
tends to zero as v/c, and that with electric field as (v/c) 2 . Moreover, the part of interaction 
(|56j) which is not related to ^-factor, is reducible in spin; in other words, the structure SiSj 
in it cannot be rewritten as an irreducible tensor s^Sj — (l/3)5.y-s 2 . The interaction (|36"D is 
not at all a quadrupole one. However, its asymptotic behaviour at 7 — > 00 is of interest. In 
this limit 



2m 2 



(g - 1) (s [v x V]) [(sB) - (sv)(vB) - (s [v x E])]. 



(37) 
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Surprisingly, the asymptotical formulae fl30| ) and (^) coincide up to a factor and a total time 
derivative. To prove this fact, it is convenient to introduce three orthogonal unit vectors 



[v x s] 

v; p = w^w\ ] c = [vxp] - 

Using the completeness of this basis and the equation E = [V x B] , and neglecting as well 
a total derivative in t , one can easily check that 

[(sV) - (vs)(vV)] [(sE) - (sv)(vE) + (s[v x B])] 

= [vxs] 2 (CV)[(CE) + (pB)], 

coincides indeed with 

(s [v x V]) [(sB) - (sv)(vB) - (s [v x E])] 

= - [v x s] 2 (pV) [(p[v x B]) + (pE)] . 
Thus, there is a special value of the bare quadrupole moment ( |2"9"D 

es 2 1 
Q = - 2 ^" 1 )^' or ^ = ^-^^2 (38) 

(let us recall that we consider now a classical situation, when s ^> 1), at which the total 
interaction, quadratic in spin, L 2s + L\ s , falls down asymptotically with energy. 

The situation resembles that which takes place for the interaction linear in spin. It is 



well-known (see, for instance, [10], [30], |3TJ] ) that there is a special value of g-factor, g = 2 , 
at which the interaction linear in spin decreases at 7 — > 00 . It follows immediately from 
formula (||) for the first-order Lagrangian. Thus, putting additionally g = 2 , we obtain 



es 2 



Q = -2 — , or / e , 2 = — . (39) 

Let us note that the choice g = 2 for the bare magnetic moment is a necessary (but 
insufficient) condition of renormalizability in quantum electrodynamics 0, [30], 31]. It is 



satisfied not only for electron, but also for the charged vector boson in the renormalizable 
electroweak theory. 

In some respect, however, the situation with the singled out values (|38|), ( |39] ) of the 
quadrupole moment differs from the situation with ^-factor. Conditions (|38|), fl3~9|) , as distinct 
from the condition g = 2 , are not universal, they take place only for large spins, s ^> 1 , 
in other words, they refer only to classical objects with internal angular momentum. In 
particular, for the charged vector boson of the renormalizable electroweak theory the bare 
quadrupole interaction is absent at all, f e>2 = . The quadrupole moment of this particle is 
(in our language) of the induced nature, it is given by formula (|3l|) at s = 1 g = 2 . 
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5 Spin precession in gravitational field 



In this section we present a simple and general derivation of the equations of the spin pre- 
cession in a gravitational field. This approach not only allows one to reproduce and easily 
generalize known results for spin effects. Pointed out here remarkable analogy between grav- 
itational and electromagnetic fields allows one to transcribe easily the results of two previous 
sections for the case of an external gravitational field. 

It follows from the angular momentum conservation in flat space-time taken together with 
the principle of equivalence that the spin 4- vector S* 1 is parallel transported along the particle 
world-line. The parallel transport of a vector along a geodesic x M (r) means that its covariant 
derivative vanishes: 

<«> 

Let us go over to the tetrad formalism natural for the description of spin. In virtue of relation 
([401) the equation for the tetrad components of spin S a = S^e* is 



DS a dS a 

Here 



Dt (h - S»e%X = rf b lbcd u d S c . (41) 



labc = e atl . u e^e u c = -^ bac (42) 



are the Ricci rotation coefficients |[32| . Certainly, the equation for the tetrad 4- velocity 
components is exactly the same: 

dn a 

-T= V ab lbc d u d u c . (43) 
dr 

The meaning of Eqs. (fUD, (|43|) is clear: the tetrad components of both vectors vary in the 
same way, due to the rotation of the local Lorentz vierbein only. 

In the exactly same way, the 4-dimensional spin and velocity of a charged particle with 
the gyromagnetic ratio g — 2 precess with the same angular velocity in an external electro- 
magnetic field, in virtue of the Bargman, Michel, Telegdi equation 0, |25| (at g = 2) and the 
Lorentz equation: 

dSa e p Q b. e p „,& 

— rr — — ) — rr — — r ah u . 
dt m at m 

Therefore, the evident correspondence takes place: 

— Fab < ► labcU C . (44) 

m 

It allows one to obtain the precession frequency u> of the 3-dimensional spin vector s in 
external gravitational field from expression @ by means of simple substitution 

6 1 6 

— Bi — ► - -e ik a Hc u c ; —Ei — ► ~/ oic u c '. (45) 
ml m 
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This frequency equals 

( I u k 

u i = ~ e iki ( iflkic + u0 - 1 7o; c ) • (46) 



u° 



The overall factor 1/it^ in this expression is related to the transition in the lhs of Eq. (41) 
to the differentiation in the world time t: 

d dt d q d 
dr dr dt w dt 

The quantity v? w is supplied with the subscript w to emphasize that this is a world component 
of 4-velocity, but not a tetrad one. All other indices in expression ( [46]) are tetrad ones, 
c = 0, 1, 2, 3; i,k,l = 1, 2, 3. The corresponding spin-dependent correction to the Lagrangian 
is 

Lug = suj . (47) 

As an illustration of formulae ( Pf ) , ( |4"7| ) , let us apply them to the cases of spin-orbit and 
spin-spin interactions. We will restrict, as it is usual in the problem discussed, to the linear 
approximation in the gravitational field. However, in our approach, as distinct from the 
standard ones, both problems admit an elementary solution at an arbitrary particle velocity. 

Tetrads e a(U are related to the metrics as follows: 

To linear approximation one can put g^ v = r]^ u + h^ v and do not distinguish between the 
tetrad and world indices in e afJj . The ambiguity in the choice of tetrads will be fixed by 
choosing the symmetric gauge e^ u = e u ^ . Then 

1 

Using expression ( f4~2"D for the Ricci coefficients, we find to linear approximation 

labc = ^(h bc ,a- h ac ,b)- (48) 

Let us start with the spin-orbit interaction. In the centrally symmetric field created by a 
mass M, the metrics is 

2kM 2kM 
hoo — — ; h mn — — o mn - (49) 

Here nonvanishing Ricci coeffeicients are 

lijk = -^-{Vjk r i ~ dikTj) , 7oj = — —n ■ (50) 

Their substitution into formula ( |46"1) gives the following expression for the precession fre- 
quency: 



27+I kM 

uj 1s = — — — — v x r. 51 

7+1 r 6 
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In the limit of low velocities, 7 — > 1, the answer goes over into the classical result 

We go over now to the spin-spin interaction. Let the spin of the central body be s . 
Linear in s components of metrics, which are responsible for the spin-spin interaction, are: 

, [so x rL 
h 0i = 2k l0 r3 U . 

Here nonvanishing Ricci coefficients equal 

, ( [s x r}j [s x r]i [s x t]j 

Hjo = fc(Vi r3 - V, r3 ) , 7 oij = -fcV, r3 . (52) 

The frequency of the spin-spin precession is 

u ss = -k (2 - - ] (s V)V- 

+ k -J— [v(s V) - s (vV) + (vs ) V] (W) - . (53) 
7+1 r 

In the limit of low velocities this formula also goes over into the corresponding classical result 

I- 

In the conclusion of this section let us note that in the case of an external gravitational 
field there is no covariant expression for the force linear in the particle spin. In other words, 
the deviation from geodesies of the trajectory of a spinning particle is not described by the 
Riemann tensor. In this case there is a unique possible covariant structure, up to a factor (in 
Ref. g] it equals -l/2m): R^ ab u u S ab . As it was mentioned already in Introduction, that 
covariant description (as distinct from our our formulae (]46|), ([47])) contradicts the classical 
results in the limit of low velocities. 



6 Equations of motion of spinning particle 
in gravitational field. General approach 

Equations of motion in an external gravitational field to any order in spin are constructed 
analogously to the equations of motion in the case of electromagnetic field. 

We write down the scattering amplitude in a weak external gravitational field h^ u as 

- \T^hT. (54) 

The matrix element T^ u of the energy-momentum tensor between the states of momenta k 
and k' can be written as 

T^ u = 4){k') l^p,jp v F i+\ (P^ux + Pu^^x) q X F 2 
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+ S v q,){Sq) + V^(Sq) 2 ] F A ) ip(k). (55) 
The scalar operators Fi in this expression are also expanded in powers of r = (Sq) 2 : 

F < (t,r) = X;/i,2n(t)r n . 

n=0 

It is easy to convince oneself that the total number of invariant form factors fi^n is 4s + 2 
and 4s + 1 for integer and half-integer spin, respectively. The independence of the four tensor 
structures entering expression ( [55] ) is obvious. As to the completeness of the expansion, it 
can be checked, for instance, by demonstrating that the calculation of the total number of 
invariant form factors in the annihilation channel leads to the same result as the above one. 

In the generally covariant form the structure {f]^q 2 — q^q u ) h^ v corresponds to the scalar 
curvature R, and [S^S u q 2 — (S^q u + S u q^)(Sq) + rj^ u (Sq) 2 } W corresponds to the product 
R^S^S" , where R^ v is the Ricci tensor. Since we are interested in the equations of motion 
in a sourceless field, the corresponding terms in the expansion (|55|) will be omitted. 

As well as in electrodynamics the charge conservation dictates the condition / e ,o(0) = 1 , 
here the energy conservation leads to /i,o(0) = 1 . As to the term in the amplitude ([54]) 
which contains /2,o > it is convenient to rewrite it otherwise, using the analogy ( |4"4| ) with 
electromagnetic field. Substituting into the corresponding electromagnetic term 

i^-^(k')i: ab F ab ^(k) 

oe 

g = 2, (e/m)F a b = f a b = lab c u c , we arrive at the following contribution to the Lagrangian of 
gravitational interaction: 

• 1 4,(k')i; ab f ab ij(ky, (56) 



4«o 

here, as usual, u° w = e/m. Using for 7 a {, c the linear approximation (j48f) , one can check easily 
that the expression (|^) corresponds indeed to the discussed contribution to the amplitude 
at the condition /2,o = 1- Thus, in gravity the value of one more form factor at zero 
momentum transfer t is fixed. It corresponds to the conservation of angular momentum. 
This circumstance was pointed out many years ago in Ref. [|33|1. 

Let us come back to the convection term in formula (|54]) . As it was the case in electro- 
dynamics, here when going over to spinors in the rest frame, the term of first order in spin 
is written as 

8e m u° + 1 
Using relations (|25"D , (|48|), we obtain 



[s[vxV])V (57) 



P^VitV -> -p^i-dkh^ + d^h kv ) - 2p a -f ak 
Thus the expression ( p?[ ) is presented through the Ricci coefficients: 

1 u° 



e mnk s m v n u a u c ^ 

vP w U {> + 1 



As can be easily checked, the expressions ([56]) and (|58|) in sum reproduce the Lagrangian 
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7 Equations of motion of spinning particle 
in gravitational field. Second order in spin 

Let us investigate now the effects of second order in spin in the equations of motion in 
a gravitational field. In the case of a binary star these effects are of the same order of 
magnitude as the spin-spin interaction at comparable spins of components of the system 



22]. The influence of the latter on the characteristics of gravitational radiation becomes 



noticeable for a system of two extreme black holes [|19|| . Correspondingly, the spin effects 
of second order in the equations of motion get appreciable if at least one component of a 
binary is an extreme black hole |2^]. Therefore, the investigation of these effects is not of 
a purely theoretical interest only. In principle they can be observed with the gravitational 
wave detectors under construction at present. 

An obvious source of the second order spin effects is the term 

L 2sg = -h,2^-p»p u {Sq) 2 h, v (59) 
oe 

in the amplitude (|54]). Due to the relation 

P ll p v q a qf3h lJLU = p fl p v {q a qph^ v + q^h^ - q a q u h^ - q^h ua ) -> 2^fR mvfi , 
the Lagrangian ( |59"D is rewritten via the Riemann tensor: 



L 2sg = ~^u a S b u c S d R ahcd . (60) 



Instead of /i 5 2 , we have introduced a dimensionless parameter k: 

K 



1.2 



2m? 



Then it is convenient to use the Petrov representation for the components of the Riemann 
tesor (see Ref. P^fl): 



Ekl — RokOl, Ekl — Elk', C k \ — ^kmn^lrsRmnrs , Cm — Clk] 

Bm = le lrs R 0krs , B kk = 0. (61) 

We restrict to the case of a sourceless gravitational field. Then, at R a b = 0, additional 
simplifications take place: 

Cm = —Em, Bm = B ik , E kk = C kk = 0. (62) 

Finally, we arrive at the following interaction Lagrangian quadratic in spin: 



K 

J 2sg = 



2e 



(2u 2 + 1) E k i - 2 ( 2 — - ) u k u m Ei m + 5 k iu m u n E„ 

u Q + 1/ 
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+ 



— 2 U ekmnU m B n l + 



K + i) 2 

2 



r U} t UiU rn 'lL n E r) 



(63) 



M + 1 

To avoid misunderstanding, we note that all three-dimensional indices in this equation (and 
next two ones) are in fact contravariant. 

As well as in electrodynamics, here along with the "bare" interaction (|63|) , there is an 
induced one quadratic in spin. Its explicit form can be obtained most easily by putting 
in electromagnetic formula fl36|) g = 2 and by making the substitution (45). We take into 
account also the correspondence 

qi7abcU C = (qaabc- qclabi)U C -> l{di^ abc - d c ^ abi ) U° -> 1 R abci U C . 

Finally, using relations flBIp, (62]), we obtain the following result for the induced interaction: 



U 



1 



E 



1 



1 



2 Uo 



K + i) 2 
i 

M + 1 



kl^m^n-Emn 



U + 1 
1 



(«0 



1 



(«o + I) 2 

UkUlUmUnE, 



U k U m E; 



I in 



(64) 



^kmnUmB 



Til 



u + 1 

As well as in the electromagnetic case, the induced interaction here tends to zero in the 
nonrelativistic limit ~ v/c, and the spin factor in it, s^si , is not an irreducible tensor. 

The asymptotical behaviour of L 2sg and L l 2sg is the same: both Lagrangians increase 
linearly with energy. However, in this case also the coefficient in the "bare" interaction can 
be chosen in such a way, k = 1 , that the total Lagrangian of second order in spin decreases 
(as well as the analogous interaction in electrodynamics) when energy tends to infinity. At 

K = 1 



+ L 



2sg 



e(u° + 1) 



u E 



1 



A7 



U 



UkU m Ei m 



+ 



2(uo + 1 



(65) 



8 Gravimagnetic moment. Multipoles of black holes 

There is a deep analogy between the linear in spin interaction of the magnetic moment with 
the electromagnetic field 

Cem = ~^-E ab S ab (66) 

Am 

and the "bare" gravitational Lagrangian ( pD| ) quadratic in spin EDI. (It is more convenient 
here to write the gravitational Lagrangian, like C em , for the proper time r , and not for 
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the world time t, i.e., to multiply expression (|60|) by e/m.) This analogy is based on the 
following observation. It is well-known that the canonical momentum id^ enters relativistic 
wave equations for a particle in electromagnetic and gravitational external fields via the 
combination 

n M = id^ - eA^ - h: ab lab ^. 

It follows from the structure of the commutator (or Poisson brackets in the classical limit) 

[n^ry = -i(eF^- h: ab R ah ^) 

that in a sense — ^T, ab R ab ^ u plays the same role in gravity as eF^ u in electromagnetism. It 
is quite natural then that the gravitational analogue of the electromagnetic spin ineraction 
©is 

£gm = -^Rabcd,S a S° . (67) 

One can easily demonstrate that expressions (^) and (^) coincide (up to a factor e/m). 
It is sufficient to this end to take into account the relation S ab = e abcd S c Ud , as well as the 
identity 

fS 1 ef gh p p 

fl-abcd ^ ^ab ^cd efgh ^abcd > 

which is valid for a sourceless gravitational field. 

It is natural to define in analogy with the magnetic moment 

2m 

the gravimagnetic moment 

— ———S ab s cd 
2m 

The gravimagnetic ratio k, like the gyromagnetic ratio g in electrodynamics may have in 
general any value. Still, it is natural that in gravity the value k = 1 is as singled out as g = 2 
in electrodynamics. In any case, at g — 2 and k = 1 the spin equations of motion have the 
simplest form. 

For a classical object the values of both parameters g and k depend in general on its 
various properties. However, for black holes the situation is different. It follows from the 
analysis of the Kerr-Newman solution that the gyromagnetic ratio of a charged rotating black 
hole is universal (and the same as that of the electron!): g = 2 f34 . 



We will demonstrate that for the Kerr black hole the gravimagnetic ratio is k = 1. This 
result follows in fact from the analysis of the motion of spin of a black hole in an external 
field carried out in Ref. [0 (though this statement by itself was not explicitly formulated 
there). We will present here an independent and, in our opinion, more simple derivation of 
this important result. 

We will consider the interaction (^) at large distance from a Kerr hole at rest. To linear 
approximation in the created gravitational field the corresponding Lagrangian density can 
be written as 



L= flsxVlJsxV]^). (6* 
4m 
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The thus induced correction to the energy-momentum tensor has space components only: 

STy = £-[a x V]i[s x V]^(r) . (69) 
Let us find the corresponding correction to the metrics. In the gauge 

h^, v = o, v= v- \vv»K ( 7 °) 

the static Einstein equations in the linear approximation are 

The induced corrections to metrics equal 

/ioo = k-(sV) 2 -, (71) 

m r 



h i:j = — k — {2 [s x V]i[s x V]j + ^(sV) 2 } - . (72) 



fc 
rn 

Let us compare now the obtained corrections (|7T]) , (|?2|) with the corresponding contribu- 
tion to the Kerr metrics. In the Boyer-Lindquist coordinates this metrics is 

ds l= (l_^) dt 2_ ^ r 2_ E ^2_ (r 2 +a 2 + !^ sin2 ^2 sin 2 e+ ^ sin 2^ t) (?3) 

Zj ZA Zj Zj 

where A = r 2 — r 9 r + a 2 , S = r 2 + a 2 cos 2 # , a = s/m. At r 3 = the metrics (|73|) 
describes a flat space in spheroidal coordinates p2| . Meanwhile, it is Cartesian coordinates 
which correspond in the flat space to the gauge ([701). The transition from the spheroidal 
coordinates to Cartesian ones is carried out with the required accuracy by the substitution 



r + 



afar) — ra 2 



2r 2 

In the Cartesian coordinates the spin-dependent part of the 00-component of the metrics 

,oo=l-^ + ^(3cos 2 *-l) 

coincides evidently with h 00 from formula (|7T| ) at k — 1. For the space components of the 
Kerr metrics as well their spin-dependent part reduces to the expression (72) at k — 1. 

Let us note that the motion of the Kerr black hole in an external gravitational field is 
not described by the Papapetrou equation even if one leaves aside the problem of spin-orbit 
interaction linear in spin. The point is that this equation refers to the case k = JT5 . 



It is proven in the same way that for a charged Kerr hole as well the gravimagnetic ratio 
n = 1. Moreover, it can be proven that the electric quadrupole moment of a charged Kerr 
hole equals also that value, 



Q = -2 e 4 
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at which the interaction quadratic in spin decreases with energy. One can demonstrate 
that other, higher multipoles of a charged Kerr hole as well possess just those values which 
guarantee for an interaction of any order in spin (but of course, linear in an external field) 
the asymptotic decrease with energy. 

We are grateful to LA. Koop, R.A. Sen'kov, and Yu.M. Shatunov for useful discussions. 
The work was supported by the Russian Foundation for Basic Research through Grant No. 
95-02-04436-. 
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